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Abstract. From the first epoch observations of the VIMOS VLT Deep Survey (VVDS) up to z — 1.5 we have 
derived luminosity functions of different spectral type galaxies. The VVDS data, covering ~ 70% of the life of the 
Universe, allow for the first time to study from the same sample and with good statistical accuracy the evolution 
of the luminosity functions by galaxy type in several rest frame bands from a purely magnitude selected sample. 
The magnitude limit of the VVDS (Iab ~ 24) is significantly fainter than the limit of other complete spectroscopic 
surveys and allows the determination of the faint end slope of the luminosity function with unprecedented accuracy. 
Galaxies have been classified in four spectral classes, from early type to irregular galaxies, using their colours and 
redshift, and luminosity functions have been derived in the U, B, V, R and I rest frame bands for each type, in 
redshift bins from z = 0.05 to z = 1.5. In all the considered rest frame bands, we find a significant steepening of 
the luminosity function going from early to late types. The characteristic luminosity M* of the Schechter function 
is significantly fainter for late type galaxies and this difference increases in the redder bands. Within each of the 
galaxy spectral types we find a brightening of M* with increasing redshift, ranging from 0.5 mag for early type 
galaxies to ~ 1 mag for the latest type galaxies, while the slope of the luminosity function of each spectral type 
is consistent with being constant with redshift. The luminosity function of early type galaxies is consistent with 
passive evolution up to z ~ 1.1, while the number of bright (Mb ab < —20) early type galaxies has decreased by 
~ 40% from z ~ 0.3 to z ~ 1.1. We also find a strong evolution in the normalization of the luminosity function 
of latest type galaxies, with an increase of more than a factor 2 from z ~ 0.3 to z ~ 1.3: the density of bright 
(Mb ab < —20) late type galaxies in the same redshift range increases of a factor ~ 6.6. These results indicate a 
strong type-dependent evolution and identifies the latest spectral types as responsible for most of the evolution 
of the UV-optical luminosity function out to z — 1.5. 
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1. Introduction 

An unbiased and detailed characterization of the luminos- 
ity function (LF) of field galaxies is a basic requirement 
in many extragalactic issues. 

At present the local luminosity function is well constrained 
by the results obtained by the 2dF Galaxy Redshift Survey 
(2dFGRS, Norberg et al. EMU and by the Sloan Digital 
Sky Survey (SDSS, Blanton et al. These surveys 

measure redshifts for 10 5 — 10 6 galaxies over a large area, 
and therefore explore well the properties of the local 
(z < 0.3) Universe. Such large numbers of objects also 
allow to study the luminosity functions (as well as the 
correlation functions and other properties) for galaxies 
of different types, defined on the basis of colours and/or 
spectral properties. A critical analysis of the luminosity 
functions depending on galaxy type as measured from the 
various redshift surveys, as well as a comparison of the 
different results, can be found in de Lapparent (2003). 
Madgwick et al. <12002)l . analyzing 2dFGRS data, find a 
systematic steepening of the faint end slope and a fainten- 
ing of M* of the luminosity function as one moves from 
passive to active star forming galaxies. Similar results are 
found by Blanton et al. I|2001[) for the SDSS sample, mov- 
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Fig. 1. Rest frame U — B (upper panel) and B — I (lower 
panel) colours for galaxies of different types: type 1 (solid 
line), type 2 (dotted line), type 3 (short dashed line), type 
4 (long dashed line). 



ing from the redder to the bluer galaxies. 
For what concerns the high redshift Universe, several 
studies in the past ten years have aimed to map the 
evolution of the luminosity function. However, because 
of the long exposure times required to obtain spec- 
tra of high redshift galaxies, spectroscopic surveys were 
limited to a few 10 2 objects. The Canadian Network 
for Observational Cosmology field galaxy redshift survey 
(CNOC-2, Lin et al. lTOOOl) and t he ESO Sculptor Survey 
(ESS, de Lapparent et al. 12003(1 derived the luminosity 
function up to z ~ 0.5 using ~ 2000 and ~ 600 redshifts, 
respectively, de Lapparent et al. I|2003|l find a behaviour 
of the LF by type similar to the local one derived from 
2dFGRS and a strong evolution of a factor 2 in the vol- 
ume density of the late type galaxies with respect to the 
early type galaxies. Lin et al. (1999) find for early type 
galaxies a positive luminosity evolution with increasing 
redshift, which is nearly compensated by a negative den- 
sity evolution. On the contrary, for late type galaxies they 
find a strong positive density evolution, with nearly no lu- 
minosity evolution. At higher redshift, the Canada France 
Redshift Survey (CFRS, Lilly et al. 11005)1 allowed to study 
the luminosity function up to z ~ 1.1 with a sample of 
~ 600 redshifts. From this survey, the LF of the red pop- 
ulation shows small changes with redshift, while the LF 
of the blue population brightens by about one magnitude 
from z ~ 0.5 to z ~ 0.75. Other results suggest a strong 
number density evolution of early type galaxies (Bell et 
al. 120041 Faber et al. l2006[> : conversely, in the K20 survey 
(Cimatti et al. 12002)) . Pozzetti et al. <j2003")l found that red 
and early type galaxies dominate the bright end of the 
LF and that their number density shows at most a small 
decrease (< 30%) up to z ~ 1 (see also Saracco et al.^006 
and Caputi et al. 2006). 

Luminosity function estimates at higher redshift and/or 
with larger samples are up to now based only on photo- 
metric redshifts, like the COMBO-17 survey (Wolf et al. 
12005)1 and the analysis of the FORS Deep Field (FDF, 
Gabasch et al. 12001)1 and the Hubble Deep Fields (HDF- 
N and HDF-S, see e.g. Sawicki et aL lTOOTI Poli et al. 12001 
2003}); most of these projects derived also the luminosity 
function for different galaxy types. Wolf et al. (P003 ) find 
that early type galaxies show a decrease of a factor ~ 10 
in 4>* up to z = 1.2. Latest type galaxies show a bright- 
ening of about one magnitude in M* and a <jf increase of 
a factor ~ 1.6 in their highest redshift bin (z ~ 1.1) in 
the blue band. Giallongo et al. 1|2005J) . using HDFs data, 
find that the B band number densities of red and blue 
galaxies have a different evolution, with a strong decrease 
of the red population at z = 2 — 3 and a corresponding in- 
crease of the blue population. Dahlen et al. : 2005;. using 
GOODS data, claim that the starburst population frac- 
tion increases with redshift by a factor of 3 at z = 2 in 
the U band. 

Although photometric redshifts represent a powerful tool 
for deep surveys, their precision strongly relies on the 
number of used photometric bands, on the templates and 
on the adopted training procedure; moreover, they are af- 
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Fig. 2. Co-added VVDS spectra of galaxies of the four different types in various redshift bins, shown in rest frame 
wavelength. The different types are indicated with different line strengths: the ligther the line, the later the galaxy 
type. The redshift bin is indicated in the label of each panel. The flux on the y-axis is in arbitrary units and the 
spectra of the various types are arbitrarily rescaled for clarity. 



fected by the problem of "catastrophic errors" , i.e. objects 
with a large difference between the spectroscopic and the 
photometric redshift. 

A major improvement in this field is obtained with the 
VIMOS VLT Deep Survey (VVDS, Le Fevre et al. l2003bl> 
and the DEEP-2 Galaxy Redshift Survey (Davis et al. 
2003). The VVDS is an ongoing program to map the evo- 
lution of galaxies, large scale structures and AGNs from 
the redshift measurements of ~ 10 5 objects down to a 



magnitude Iab — 24, in combination with a multiwave- 
length dataset from radio to X-rays. 

From the analysis of the evolution of the global luminos- 
ity function from the first epoch VVDS data (Ilbert et 
al. I2005f) . we found a significant brightening of the M* 
parameter in the U, B, V, R and I rest frame bands, go- 
ing from z — 0.05 to z = 2. Moreover, we measured an 
increase of the comoving density of bright galaxies: this 
increase depends on the rest frame band, being higher in 
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Fig. 3. Observed fraction of bright galaxies (Mb ab — 
5log(h) < —20) of different types as a function of redshift. 
Error bars are la Poisson errors. 



the bluest bands. 

Among the other results of this survey, we recall the study 
of the radio selected objects (Bondi et al. I2003[) and of 
their optical counterparts (Ciliegi et al. I2005f) . the evolu- 
tion of the clustering properties (Le Fevre et al. 2005b 
Polio et al. 12.005) and of the bias parameter (Marinoni 
et al. EDIEl). Moreover, from the joined GALEX-VVDS 
sample, we derived the evolution of the far UV luminos- 
ity function (Arnouts et al. 2005) and luminosity density 
(Schiminovich et al. l2005[) . 

In this paper we study the evolution of the luminosity 
functions of galaxies of different spectral types based on 
the VVDS data. This sample allows to perform this anal- 
ysis for the first time with excellent statistical accuracy 
over a large redshift range (0.05 < z < 1.5). 
The plan of the paper is the following: in sect. 2 we briefly 
present the first epoch VVDS sample, in sect. 3 we de- 
scribe the galaxy classification and in sect. 4 we illustrate 
the method we used to estimate the luminosity functions. 
In sect. 5 we compare the luminosity functions of the dif- 
ferent galaxy types and in sect. 6 we show the evolution 
with redshift of the luminosity functions by type. Finally 
in sect. 7 we compare our results with previous literature 
estimates and in sect. 8 we summarize our results. 
Throughout the paper we adopt the cosmology Q m = 
0.3 and fl A = 0.7, with h = H /100 km s" 1 Mpc" 1 , 
Magnitudes are given in the AB system and are expressed 
in the five standard bands U (Bessel), B and V (Johnson), 
R and I (Cousins). 



2. The first epoch VVDS sample 

The VVDS is described in detail in Le Fevre et al. 
(2005a): here we report only the main characteristics of 
the sample used for the analysis presented in this paper. 
The entire VVDS is formed by a wide part on 4 fields 
(which is not used in this paper), and by a deep part, 
with spectroscopy in the range 17.5 < Iab < 24 on the 
field 0226-04. Multicolour photometry is available for 
each field (Le Fevre et al. 12004a}) : in particular, the B, V, 
R, I photometry for the 0226-04 deep field is described 
in detail in McCracken et al. (2003). Moreover, U band 
(Radovich et al. I2004J) and J and K band (Iovino et al. 
2005) data are available for smaller areas of these fields. 
Starting from these photometric catalogues, spectroscopic 
observations were performed with the Visible Multi- 
Object Spectrograph (VIMOS, Le Fevre et al. I2003a|l 
mounted on the ESO Very Large Telescope (UT3). The 
selection of objects for spectroscopic observations was 
based only on magnitude, without any other colour or 
shape criteria. 

Deep spectroscopic observations (17.5 < Iab < 24) 
were performed also on the Chandra Deep Field South 
(VVDS-CDFS, Le Fevre et al. l2004hll . starting from the 
EIS I band photometry and astrometry (Arnouts et al. 
IMTTjl . Multicolour U, B, V, R and I photometry for this 
sample is available from the COMBO- 17 survey (Wolf et 

ai. riom 

Spectroscopic data were reduced with the VIMOS 
Interactive Pipeline Graphical Interface (VIPGI, 
Scodeggio et al. 12717151 Zanichelli et al. l2"TOHjl and 
redshift measurements were performed with the KBRED 
package (Scaramella et al. 12006|l and then visually 
checked. Each redshift measurement was assigned a 
quality flag, ranging from (failed measurement) to 4 
(100% confidence level); flag 9 indicates spectra with 
a single emission line, for which multiple solutions are 
possible. Further details on the quality flags are given in 
Le Fevre et al. H2005afl . 

The analysis presented in this paper is based on the 
first epoch VVDS deep sample, which has been obtained 
from the first observations (fall 2002) on the fields 
VVDS-02h and VVDS-CDFS, which cover 1750 and 450 
arcmin 2 , respectively. We eliminated from the sample 
spectroscopically confirmed stars and broad line AGNs, 
remaining with 6477 4- 1236 galaxy spectra with secure 
spectroscopic identification (flag 2, 3, 4, 9), corresponding 
to a confidence level higher than 75%. Redshifts with 
flags and 1 are taken into account statistically (see sect. 
4). This spectroscopic sample, which is purely magnitude 
selected, has a median redshift of ~ 0.76. 

3. Galaxy classification 

Galaxies have been classified using all the multicolour in- 
formation available; in the VVDS-02h field B, V, R and 
I band magnitudes are available for all galaxies, while U 
band data are available for 83% of the galaxies. For the 
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VVDS-CDFS sample U, B, V, R and I photometry from 
the COMBO-17 survey is used. 

Absolute magnitudes are computed following the method 
described in the Appendix of Ilbert et al. {2005). The K- 
correction is computed using a set of templates and all the 
photometric information (UBVRI) available. However, in 
order to reduce the template dependency, the rest frame 
absolute magnitude in each band is derived using the ap- 
parent magnitude from the closest observed band, red- 
shifted at the redshift of the galaxy. With this method, 
the applied K-correction is as small as possible as possi- 
ble. 

For each galaxy the rest frame magnitudes were matched 
with the empirical set of SEDs described in Arnouts et 
al. I|1999[) . composed of four observed spectra (CWW, 
Coleman et al. ll98Uf> and two starburst SEDs computed 
with GISSEL (Bruzual & Chariot ITMS)l . The match is 
performed minimizing a \ 2 variable on these templates at 
the spectroscopic redshift of each galaxy. The same proce- 
dure has been applied by Lin et al. (jl999|) to the CNOC-2 
survey up to z ~ 0.55. This approach is also similar to 
that adopted by Wolf et al. PJ03|| for the COMBO-17 
survey, but we have the advantage of using spectroscopic 
redshifts, while they had to rely on photometric redshifts. 
Galaxies have been divided in four types, corresponding 
to the E/SO template (type 1), early spiral template (type 
2), late spiral template (type 3) and irregular template 
(type 4). These types are based on the four CWW tem- 
plates: type 4 includes also the two starburst templates. 
The numbers of galaxies for each type are listed in Table 
□ 

In order to have an idea of the correspondence of these 
types with colours, we report here the rest frame colours 
for each template: type 1, 2, 3 and 4 have Bab — 
Iab =1-58, 1.11, 0.79 and 0.57, respectively. Given these 
colours, a rough colour subdivision for each class is 1.3 < 
Bab - Iab for type 1, 0.95 < B A b - Iab < 1-3 for type 2, 
0.68 < B AB -Iab < 0.95 for type 3 and B A b~Iab < 0.68 
for type 4. However, we remind that these colour ranges 
are only indicative and our classification scheme is based 
on the whole multicolour coverage. In Fig^we show the 
U — B and B — I colour distributions for the galaxies of 
our sample divided according to type. From this figure 
it is clear that, although the different types have differ- 
ent colour distributions, they present significant overlaps. 
This fact is a consequence of classification schemes using 
template fitting on multicolour data. 
Note that, in order to avoid to be model dependent, we 
did not apply to the templates any correction aimed at 
taking into account colour evolution with redshift. It is 
well known that the colour of a simple stellar population 
subject to passive evolution was bluer in the past. In prin- 
ciple, this could imply that galaxies classified as type 1 at 
low redshift might be classified differently at higher red- 
shift. Indeed, this effect has been invoked by the authors 
who found negative evolution in the luminosity function 
of "red" galaxies (see f.i. Wolf et al. 2003). In order to ver- 
ify this hypothesis, we applied our classification scheme to 



Table 1. Numbers of galaxies of different types in the 
sample 



Type 




total 


7713 


type 1 


730 


type 2 


1290 


type 3 


2622 


type 4 


3071 



synthetic spectra (Bruzual & Chariot 1993) of ellipticals 
(i.e. simple stellar populations and exponentially declining 
star formation with time scales of 0.1 Gyr and 0.3 Gyr) 
with formation redshift between Zf orm = 2 and 20. We find 
that all ellipticals with Zf orrn > 2 would be classified as 
type 1 objects, even at z ~ 1. 

In order to check, at least on a statistical basis, the consis- 
tency between this photometric classification and average 
spectral properties, we summed the normalized spectra of 
all galaxies in each of the four types. The resulting average 
spectra are shown in Fig|3for each type in various redshift 
bins. This figure confirms the robustness of our classifica- 
tion scheme: moving from type 1 to type 4 objects, the 
composite spectra show an increasingly blue continuum, 
with emission lines of increasing strength. This confirms 
that the four types show different spectral features and 
therefore represent different classes of objects. 
For the VVDS-CDFS sample, HST-ACS images are avail- 
able. Using these data, Lauger et al. (|2006ll classified the 
galaxies in this sample using an asymmetry-concentration 
diagram. Plotting our type 1 galaxies in this diagram, we 
find that ~ 91% of them lie in the region of bulge dom- 
inated objects, showing an excellent consistency also be- 
tween our photometric classification and a morphological 
one. 

In FigOl we plot the observed fraction of bright galaxies 
of each type as a function of redshift. We selected objects 
with Mb ab — 5log(h) < —20 because these galaxies are 
visible in the whole redshift range. From this figure it is 
clear the growing importance of bright late type objects 
with increasing redshift and the corresponding strong de- 
crease of the fraction of bright early type galaxies. 

4. Luminosity function estimate 

Luminosity functions were derived using the Algorithm for 
Luminosity Function (ALF), a dedicated tool which uses 
various estimators: the non-parametric 1/V max (Schmidt 
GUSH, C+ (Lynden Bell [MP , SWML (Efstathiou et al. 
1988J and the parametric STY (Sandage, Tammann & 
Yahil nTTTHjl . for which we assumed a Schechter function 
fSchechter I1976fl . The tool and these estimators, as well 
as their specific use in the context of the VVDS, are de- 
scribed in detail in Ilbert et al. ( 2005). 
Ilbert et al. ( 2001J have shown that the estimate of the 
global luminosity function can be biased, mainly in its 
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Fig. 4. Luminosity functions by different types in the redshift range [0.4 — 0.9] in various rest frame bands: upper 
panels U (left) and B (right); lower panels R (left) and I (right). Note that the ranges of magnitude are different in 
the various panels. The lines represent the STY estimates for type 1 (solid), type 2 (dotted), type 3 (short dashed) 
and type 4 (long dashed) galaxies. The vertical dashed line represents the faint absolute magnitude limit considered 
in the STY estimate (see text). The shaded regions represent the 68% uncertainties on the parameters a and M* , 
whose confidence ellipses are reported in the right panels. The ellipse contours are at 68% and 90% confidence level 
(solid and dotted line respectively). The points inside the ellipses represent the best fit values for type 1 (filled circle), 
type 2 (open circle), type 3 (filled triangle) and type 4 (open triangle) galaxies. 



faint end, when the band in which it is measured is far 
from the rest frame band in which galaxies are selected. 
This is due to the fact that, because of the K-corrections, 
different galaxy types are visible in different absolute mag- 
nitude ranges at a given redshift. When computing the 
global luminosity functions (Ilbert et al. l2005f) . we avoided 
this bias by using for the STY estimate, in each redshift 
range, only galaxies within the absolute magnitude range 
where all the SEDs are potentially observable. 



Even if this bias is much less important when estimating 
the luminosity function of galaxies divided by types, we 
have, however, taken it into account. The absolute magni- 
tude limits for the STY estimate are indicated with ver- 
tical dashed lines in the figures, and in the tables where 
the best fit parameters are reported (Table |21 and we 
give both the total number of objects and the number of 
galaxies within this magnitude limit. 
In order to take into account the unknown redshifts (not 
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(a) Number of galaxies in the redshift bin 

(b) Number of galaxies brighter than the bias limit (sample used for STY estimate; see the text for details) 



observed objects and failed spectra), a weight was applied 
to each galaxy following the procedure described in detail 
in llbert et al. lj!2UU5|> . 

This weight is a combination of two different contribu- 
tions: the target sampling rate and the spectroscopic suc- 
cess rate. The target sampling rate, i.e. the fraction of 
observed galaxies, corrects for the selection effects due to 
the procedure used for the mask preparation (Bottini et 
al.[2005): to maximize the number of slits, the procedure 
tends to select objects with smaller angular size on the x- 
axis of the image, corresponding to the direction in which 
the slits are placed. As a consequence, the final spectro- 
scopic sample has a bias against large objects, which pro- 
duces a mild dependence of the target sampling rate on the 
apparent magnitude. The target sampling rate is ~ 25% 
for most of the sample and is computed as a function of the 
object size (see llbert et al. 120051 for further details). The 
spectroscopic success rate takes into account the fraction 
of objects without a good redshift determination (i.e. flags 
and 1). As shown in llbert et al. (2005), these objects are 
expected to have a different redshift distribution with re- 
spect to that of the sample with measured redshift, as con- 
firmed by the use of their photometric redshifts. Given the 
fact that the spectroscopic success rate decreases for faint 
apparent magnitudes, we derived it in four magnitude bins 
as a function of redshift (using photometric redshifts, see 
Fig. 3 in llbert et al. l2005[) . The shape of the spectroscopic 
success rate is similar in all magnitude bins, showing a 
maximum at z ~ 0.7 and two minima for z < 0.5 and 
z > 1.5. Since the number of galaxies for each type is not 



large enough to reliably estimate the spectroscopic success 
rate as a function also of the galaxy type, we have used 
the global spectroscopic success rate for all galaxy types. 
In order to check the effect of the "cosmic variance", 
i.e. variations in the luminosity function due to fluctu- 
ations in the large scale structure, we applied the fol- 
lowing test on the VVDS-02h deep area. For this field 
we derived photometric redshifts (llbert et al. I2006f) 
based on both VVDS photometry (BVRIJK) and on 
new CFHT Legacy Survey photometry (ugriz), which 
has now become available in a field covering 1 sq. deg. 
( http : / / www.cfht . hawaii . edu /Science / CFHLS / ) which in- 
cludes the 1700 arcmin 2 area covered by the VVDS spec- 
troscopic survey. Then we divided the field in two non 
overlapping regions (the sub-area where spectroscopic 
data are available and the remaining area) and we com- 
pared the luminosity distributions of galaxies in these two 
samples (^0.5 sq.deg. each), in the same redshift bins 
in which the luminosity functions were derived. In each 
redshift bin the two distributions show average differ- 
ences of the order of 10%, with some larger fluctuations 
due to Poisson statistics, without any systematic trend. 
Therefore the influence of the "cosmic variance" is ex- 
pected to be limited. 

5. Comparison of the luminosity functions of 
different types 

As a first step, we compare the luminosity functions for 
galaxies of different types, in order to see which is the 
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relative behaviour of the various populations. To perform 
this comparison, we selected galaxies in the redshift range 
[0.4 — 0.9]. About 50% of the objects of our sample are 
included in this redshift interval, covering a wide range of 
luminosities (i.e. absolute magnitudes in the B band are in 
the range [—23.7; —16.8] — 5 log(/i)). Moreover, the spec- 
troscopic success rate of our survey reaches a maximum 
at z ~ 0.7 and therefore the possible dependency of the 
estimated LF on the weighting scheme described above is 
minimized in this redshift interval. 

In FigQJwe report the luminosity functions estimated with 
the STY method in the rest frame bands U, B, R and 
I, with the corresponding confidence ellipses for the a 
and M* parameters. The values of the parameters with 
their ler errors, as well as the number of galaxies for each 
type, are reported in Table El The <f>* parameters listed 
in this table are derived adopting the density estimator of 
Efstathiou et al. ( 1988), following the procedure described 
in the Appendix of Ilbert et al. ( 2005). 
TableEJshows that our estimates are based on several hun- 
dreds of galaxies for each type, and are therefore well sta- 
tistically constrained, as can be seen also from the sizes of 
the confidence ellipses in the figures. The first, very clear 
result which appears from FigE] is the significant strong 
steepening of the luminosity functions going from early to 
late types. In all bands the power law slope steepens by 
Aa ~ 1.3 — 1.5 going from type 1 to type 4 galaxies and 
galaxies of late types are the dominant population at faint 
magnitudes. 

Systematic trends are also seen in the M* parameter. In 
the reddest bands (lower panels in Fig0J, M* is signifi- 
cantly fainter for late type galaxies and this faintening is 
particularly apparent for type 4 objects. The brighter M* 
for early type galaxies reflects the fact that most of the 
more massive objects belong to this population. 
The difference of the M* values for different types de- 
creases in the B band and disappears or even changes sign 
in the U band. This behaviour is explained by the fact that 
the luminosity in the bluer bands is dominated by the light 
of young stars, produced during the star formation activ- 
ity. Galaxies of later types, which are still actively forming 
stars, are therefore more luminous in the bluer bands. 
These results are qualitatively in agreement with previous 
results from the literature, most of which at lower redshift 
(see de Lapparent 120031 for a review of the results from 
a number of surveys in the redshift range < z < 0.6). 
In particular, in almost all surveys the luminosity func- 
tion of late type galaxies is steeper and with a fainter M* 
with respect to that of early type galaxies. However, a 
quantitative comparison with previous results is difficult, 
because of the different classification schemes adopted in 
the various surveys, the different redshift ranges and se- 
lection criteria. 



6. Evolution with redshift of the luminosity 
functions by type 

We derived luminosity functions for each type in redshift 
bins in the U, B, V, R and I rest frame bands. Given our 
multicolour coverage and the explored redshift range, the 
estimate of the absolute magnitudes in the U and B rest 
frame bands are those which require less extrapolations 
(see Appendix A and Figure A.l in Ilbert et al. 2005). 
Therefore, to limit the number of figures in the paper, we 
show the results in the B rest frame band. 
Figures El El El an( i El show the luminosity function for type 
1, 2, 3 and 4 galaxies in redshift bins, obtained with C + 
and STY methods. The luminosity functions derived with 
the other two methods {1/V max and SWML) are consis- 
tent with those shown in the figures, but are not drawn for 
clarity. The dotted line in each panel represents the fit de- 
rived in the redshift range [0.4—0.9] (see previous section), 
while the dashed line is the estimate derived by fixing the 
slope a to the value obtained in the range [0.4 — 0.9]. 
In TableEHwe report the Schechter parameters, with their 
la errors, estimated for the various redshift bins, from 
z = 0.2 to z = 1.5 for each type; as a reference, in the 
last line we give the parameters derived in the redshift 
bin [0.4-0.9]. 

We do not show the results for bins where the number 
of objects is too small (less than ~ 30) to constrain the 
parameters of the luminosity function (these are the bin 
[0.05 - 0.2] for type 1 and 2 and the bin [1.2 - 1.5] for 
type 1). Given the bright magnitude limit of the survey 
(Lib > 17.5) and the small sampled volume, bright galax- 
ies are not sampled in the redshift bin [0.05 — 0.2] and 
therefore we can not constrain the M* parameter even for 
type 3 and 4 galaxies, where the number of objects is rel- 
atively high (~ 80). For this reason we show in the figures 
the luminosity function estimates in this bin, but we do 
not report the STY parameters for this redshift range in 
Table El 

In FigEl we show the confidence ellipses of the parame- 
ters a and M* in different redshift bins for the different 
types. From this figure it is possible to see that, within 
each type, the estimated slopes a in the various redshift 
bins are always consistent (within 90% confidence level) 
with each other and with the value derived in the redshift 
range [0.4 — 0.9]. Therefore, there is no evidence of a signif- 
icant change with redshift of the luminosity function slope 
within each galaxy type. Note also that the uncertainties 
on the slope estimates become quite large for z > 1: this is 
due to the fact that, even with the faint limit (Lab < 24) 
of this survey, the number of galaxies fainter than M* is 
too low to well constrain the slope. 

In each panel of Figures El 13 an d El we draw, as a ref- 
erence, the luminosity function derived in the redshift bin 
[0.4 — 0.9] (dotted line). Comparing this curve with the 
estimates of the luminosity function in the different red- 
shift bins an evolution can be seen, strongly depending 
on the galaxy type. This evolution is particularly evident 
for type 4 galaxies: going from low to high redshift there 
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is an almost continuous brightening of M* and at fixed 
luminosity the density of these galaxies was much higher 
in the past. 

The observed evolution of the luminosity function could 
be due to an evolution in luminosity and/or in density or 
both. One of the advantages of the STY method is that 
of allowing to derive the a and M* parameters indepen- 
dently from </>*, which is not possible when one directly 
fits a Schechter function on the l/V max points. 
Given the fact that we found that a is consistent with be- 
ing constant for each type, we can fix it at the reference 
value derived in the redshift range [0.4 — 0.9] and then 
study the variations of the parameters M* and <f>* as a 
function of the redshift (see upper and middle panels of 
Fig II 0(1 . These estimates are reported in Table 01 



From Fig|n3 we can see a mild evolution of M* from the 
lowest to the highest redshift bin for each type. In par- 
ticular, this brightening ranges from ^ 0.5 mag for early 
type galaxies to ~ 1 mag for the latest type galaxies. The 
only exception with respect to the general trend is that 
of type 3 objects in the bin [0.2 — 0.4], for which the best 
fit value of the M* parameter is significantly fainter than 
expexted. The reason for this discontinuity in M* for type 

3 galaxies at low redshift is not clear. 

On the contrary, the <jf parameter shows a very different 
behaviour for type 1, 2 and 3 galaxies with respect to type 

4 galaxies. The first three types show a rapid decrease of 
<fi* at low redshifts (between z ~ 0.3 and z ~ 0.5), then <j>* 
remains roughly constant up to z ~ 0.9 and finally slowly 
decreases up to z = 1.5. 
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Fig. 6. Evolution of the luminosity function in the B-band for type 2 galaxies. The meaning of the lines and the 
symbols is the same as in Fig0 



Type 4 objects, on the contrary, show an increase in <fi* at 
low redshift, then <j>* is nearly constant up to z ~ 0.8 and 
shows a rapid increase of a factor ~ 2 at z = 1.1. Then 
there seems to be a decrease from z = 1.1 to z = 1.3. 
However, this decrease can likely be a spurious effect, 
due to the fact that in this bin the estimated M* is 
very close to the bias limit (see Sect. 4). If, for exam- 
ple, we fix M* to the value obtained in the previous red- 
shift bin [1.0 — 1.2], we derive a significantly higher value 
for 4>*i i-e. (f>* — 5.83 x lQ~ 3 h 3 Mpc~ 3 (corresponding to 
4>*/<f>* e f = 1.31). Therefore the last 4>* value for type 4 



galaxies is likely to be a lower limit of the true density 
value. 

This analysis of the different trends with redshift of the <fi* 
parameter indicates that the importance of type 4 galaxies 
is increasing with redshift. However, since M* is chang- 
ing with redshift (see above), the trends of <fi* can not 
be immediately interpreted in terms of density at a given 
absolute magnitude. For this reason, we have computed 
the density of bright galaxies as a function of redshift. 
We integrated the best fit luminosity function down to 
Mb ab — 5log(h) < —20. This limit approximately corre- 
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Fig. 7. Evolution of the luminosity function in the B-band for type 3 galaxies. The meaning of the lines and the 
symbols is the same as in Fig0 



sponds to the faintest galaxies which are visible in the 
whole redshift range. In the lowest panel of Fig^|we plot 
the density of bright galaxies of each type as a function of 
redshift. 

The main results shown in this plot can be summarized in 
the following way: 

a. the density of bright early type galaxies (type 1) de- 
creases with increasing redshift, however this decrease is 
rather modest, being of the order of ~ 40% from z ~ 0.3 
to z ~ 1.1; 

b. the density of bright late type galaxies (type 4) is 



instead significantly increasing, by a factor ~ 6.6 from 
z ~ 0.3 to z ~ 1.3. 

The behaviour of type 4 galaxies is also responsible of the 
evolution of the global luminosity function measured by 
Ilbert et al. p¥)5)) . In fact, the increasing number of both 
faint and bright type 4 galaxies leads to the steepening of 
the global LF (due to the very steep slope of type 4 LF) 
and to the brightening of M* (due to the increasing frac- 
tion of bright blue objects). This fact has been directly 
checked summing the LF of all types and comparing the 
result with the global LF estimate. 
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7. Comparison with previous literature results 

Although various estimates of the luminosity function by 
galaxy type are available in the literature, a quantita- 
tive comparison of our results with previous analyses is 
not straightforward, because of the different classification 
schemes adopted in the various surveys, the different num- 
bers of galaxy types, the different redshift ranges and se- 
lection criteria. 

The strong density evolution of late type galaxies we find 
in the redshift range [0.2 — 1.5] extends to significantly 



higher redshift the results found by de Lapparent et al. 
(2004) for their latest type at z ~ 0.5. 
Among the other surveys, CNOC-2 (Lin et al. I1999f) and 
COMBO-17 (Wolf et al. ,20031) adopted a classification 
scheme somewhat similar to ours. Lin et al. (fTM5)) di- 
vided their sample of galaxies with z < 0.55 in three 
classes (early, intermediate and late type) using CWW 
templates. For early type galaxies they found a positive 
luminosity evolution, which is nearly compensated by a 
negative density evolution. On the contrary, for late type 
galaxies they found a strong positive density evolution, 
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Table 3. STY parameters (with la errors) for different galaxy types in the rest frame B band 
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(a) Number of galaxies in the redshift bin 

(b) Number of galaxies brighter than the bias limit (sample used for STY estimate; see the text for details) 



with nearly no luminosity evolution. We see at higher red- 
shift the same trend in density, but our evolution in M* 
is contained within one magnitude for each type. Wolf et 
al. 1)2003[1 used a sample of ~ 25000 galaxies with photo- 
metric redshifts, applying a classification scheme in four 
classes similar to ours but using the Kinney et al. (1996) 
templates instead of the CWW templates. 
In Figure^Jwe compare our results (solid lines) with those 
from COMBO-17 (dashed lines). Note that, since our data 
extend always to fainter absolute magnitudes, in partic- 
ular for type 1 galaxies, our estimates of the faint end 
slope are likely to be better determined: these estimates 
are derived in each redshift bin, while the COMBO-17 
slopes are fixed to the value determined in the redshift 
range [0.2 - 0.4]. This figure shows that there are signif- 
icant differences in both shapes and evolution of the LF 
estimates. The slope of the COMBO-17 LF is flatter than 
ours for type 1 galaxies, while it is steeper for types 2 
and 3 (at least up to z = 1.0). The most significant dif- 
ference between the two surveys regards the evolution of 
type 1 galaxies, for which we do not have any evidence of 
the very strong density decrease with increasing redshift 



present in COMBO-17 data. The reason for this differ- 
ence is unclear: it could be due to the use of different 
templates in the definition of the galaxy types or to a 
degeneracy between photometric redshift and classifica- 
tion, which might affect the COMBO-17 data. Bell et al. 
(2004) explained the strong negative evolution of type 1 
galaxies as a consequence of the blueing with increasing 
redshift of elliptical galaxies with respect to the template 
used for classification. In this way, at increasing redshift 
an increasing number of "ellipticals" would be assigned a 
later type, therefore producing the detected density de- 
crease observed in COMBO-17 data. However, as already 
mentioned in sect. 3, we verified that this effect does not 
affect our classification scheme, at least up to z ~ 1 and 
for simple stellar populations with Zf orm > 2. 
As a test of this possible effect, we compared the LF ob- 
tained by adding together type 1 and 2 objects. In this 
way we can check whether the differences between VVDS 
and COMBO-17 LF of type 1 galaxies is due to the fact 
that a significant fraction of high redshift type 1 galaxies in 
COMBO-17 are classified as type 2 galaxies. This compar- 
ison is shown in Fig(EI The discrepancy between VVDS 
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Fig. 9. Confidence ellipses at 90% confidence level for the a and M* parameters of the luminosity functions reported 
in Fig[3|Sl[Z|and|Sl Different line types and weights refer to different redshift bins, described in the labels; the points 
indicate the best fit values in the redshift bins [0.2 — 0.4] (filled squares), [0.4 — 0.6] (open squares), [0.6 — 0.8] (filled 
triangles), [0.8 — 1.0] (open triangles), [1.0 — 1.2] (filled circles), [1.2 — 1.5] (open circles). The large open star indicates 
the reference value obtained in the redshift bin [0.4 — 0.9]. 



and COMBO-17 LF is now reduced, but there are still 
significant differences both in slope (being the COMBO- 
17 LF steeper than that of VVDS) and in normalization, 
especially in the highest redshift bin, where the VVDS LF 
is more than a factor 2 higher than the COMBO-17 LF. 



8. Conclusions 

In this paper we studied the evolution of the luminosity 
function of different galaxy types up to z — 1.5, using 7713 
spectra with 17.5 < Yab < 24 from the first epoch VVDS 
deep sample. 

The VVDS data allow for the first time to study with ex- 
cellent statistical accuracy the evolution of the luminosity 




Fig. 10. Evolution of the parameters <f)* (upper panel) and M* (middle panel) as a function of the redshift, for 
different galaxy types. In the lowest panel the density of bright [Mb ab — 5log(h) < —20) galaxies of different types is 
shown. The slope a is fixed to the value derived in the redshift range [0.4 — 0.9]; the suffix 'ref indicates the parameters 
estimated in the redshift range [0.4 — 0.9]. Error bars are at la. 



functions by galaxy type from relatively low redshift up to 
z = 1.5 from a single purely magnitude selected spectro- 
scopic sample. The faint limiting magnitude of the VVDS 
sample allows to measure the slope of the faint end of 
the luminosity function with unprecedented accuracy up 
toz~ 1.2. The use of spectroscopic redshifts implies low 
"catastrophic" failure rate compared to the photometric 
redshifts and therefore rare populations are sampled with 
a better accuracy, like in the bright end on the luminos- 



ity function. Moreover, the use of spectroscopic redshifts 
allow to classify galaxies avoiding a possible degeneracy 
between photometric redshift and classification. 
VVDS galaxies were classified in four spectral classes us- 
ing their colours and redshift, from early type to irregular 
galaxies, and luminosity functions were derived for each 
type in redshift bins, from z = 0.05 to z = 1.5, in the U, 
B, V, R and I rest frame bands. 

We find a significant strong steepening of the luminos- 
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Fig. 11. Comparison between VVDS and COMBO-17 luminosity functions, in various redshift bins (indicated in the 
label in each panel) and for various types. Upper panels: type 1 (left) and type 2 (right) galaxies. Lower panels: type 
3 (left) and type 4 (right) galaxies. Solid lines: VVDS estimate. Dotted lines: VVDS estimate in the redshift range 
[0.4 - 0.9], plotted as a reference. Dashed line: COMBO-17 estimate from Wolf et al. 1)2003(1 . 



ity function going from early to late types: in all bands 
the power law slope steepens by Acv ~ 1.3—1.5 going 
from type 1 to type 4. Moreover, the M* parameter of 
the Schechter function is significantly fainter for late type 
galaxies. As expected, this difference increases in the red- 
der bands, reaching ~ 1.4 mag in the I band. 
Studying the variations with redshift of the luminosity 
function for each type, we find that there is no evidence 
of a significant change of the slope, while we find a bright- 
ening of M* with increasing redshift, ranging from ^ 0.5 
mag for early type galaxies to ~ l mag for the latest type 
galaxies. We also find a strong evolution in the normal- 
ization of the luminosity function of latest type galaxies, 



with an increase of more than a factor 2 in the <p* param- 
eter going from z ~ 0.3 to z ~ 1.3. The density of bright 
(M b ab —5log(h) < —20) galaxies shows a modest decrease 
(~ 40%) for early type objects from z ~ 0.3 to z ~ 1.1; 
on the contrary, the number of bright late type galaxies 
increases of a factor ~ 6.6 from z ~ 0.3 to z ~ 1.3. 
Our results indicate that the importance of type 4 galaxies 
is increasing with redshift, with an important contribu- 
tion of both bright and faint blue objects. This fact is also 
largely responsible of the evolution of the global luminos- 
ity function measured by Ilbert et al. I|2005(l . which shows 
a brightening of M * and a steepening of a with increas- 
ing redshift. Moreover, the increasing contribution of blue 
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Fig. 12. Comparison between VVDS and COMBO- 17 lu- 
minosity functions, in various redshift bins (indicated in 
the label in each panel) for the sum of types 1 and 2. The 
meaning of the symbols is the same as in Fiallll 

galaxies has been seen in the evolution of the GALEX- 
VVDS luminosity function at 1500A(Arnouts et al. l2Tjtj5|) . 
We are therefore pinpointing that the galaxies responsible 
for most of the evolution quoted in the literature belong 
to the population of the latest spectral type. The epoch at 
which a transition between a Universe dominated by late 
type galaxies and a Universe dominated by old massive 
objects occurs is at a redshift of z ~ 0.7 — 0.8. 
The fact that type 1 galaxies show only a mild evolution 
both in luminosity (positive) and in density (negative) is 
consistent with the fact that most of the objects in this 
class are old (zf orm > 2, see sect. 3) galaxies, experienc- 
ing only a passive evolution in the explored redshift range. 
More intriguing is the density evolution of type 4 galax- 
ies, which corresponds to an increasing number of bright 
star forming galaxies towards high redshift which could be 
connected to various populations of high redshift objects 
seen in multiwavelength surveys. 
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